Multi-scale magnetic study on Ni(lll) and graphene on Ni(lll) 
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We have investigated the magnetism of the bare and graphene-covered (111) surface of a Ni single 
crystal employing three different magnetic imaging techniques and ab initio calculations, covering 
length scales from the nanometer regime up to several millimeters. With low temperature spin- 
polarized scanning tunneling microscopy (SP-STM) we find domain walls with widths of 60 - 90 nm, 
which can be moved by small perpendicular magnetic fields. Spin contrast is also achieved on the 
graphene-covered surface, which means that the electron density in the vacuum above graphene 
is substantially spin-polarized. In accordance with our ab initio calculations we find an enhanced 
atomic corrugation with respect to the bare surface, due to the presence of the carbon orbitals 
and as a result of the quenching of Ni surface states. The latter also leads to an inversion of spin- 
polarization with respect to the pristine surface. Room temperature Kerr microscopy shows a stripe 
like domain pattern with stripe widths of 3 - 6 ^m. Applying in-plane-fields, domain walls start 
to move at about 13 mT and a single domain state is achieved at 140 mT. Via scanning electron 
microscopy with polarization analysis (SEMPA) a second type of modulation within the stripes is 
found and identified as 330 nm wide V-lines. Qualitatively, the observed surface domain pattern 
originates from bulk domains and their quasi-domain branching is driven by stray field reduction. 

PACS numbers: 75.70.Rf, 07.79. Cz, 61.48.Gh, 75.60.Ch, 73.20.Hb, 73.40.Ns, 81.05.ue. 



I. INTRODUCTION 

The knowledge of the magnetic properties of 3d ma- 
terials is crucial for the design of micromagnetic devices 
and the tailoring of their properties. In particular, the 
interest in Ni(lll) was recently renewed from both the 
experimental and theoretical sides. The Shockley type 
spin-split surface state of Ni(lll) was reported to play 
an important role for the magnetic properties of the 
surface However, even though bulk Ni(lll) is a con- 
ventional ferromagnet, its domain structure is not well 
investigated. Since there are no easy magnetization axes 
in the plane, a complex magnetic pattern is expectedi^ 
A model of a multiple quasi-domain branching was pro- 
posed for a Ni(lll) platelet by Hubert and Schaferi^ For 
a long time, only bitter technique data has been avail- 
able, which gives only a rough idea of the surface domain 
structure»i Ni(lll) was more recently studied by Mag- 
netic Force Microscopy (MFM), where domains of the 
order of 500 nm were observed at T = 8 K and domain 
walls shifted in perpendicular fields of i? = 25 mT4 Still, 
a complete picture of the microscopic domain pattern is 
missing. 

Ni(lll) has become again the focus of current research 
due to its role as a perfect substrate for the growth of 
graphenei^ Owing to the very small lattice mismatch, 
it grows pseudomorphically, and it has been recently 
shown that, due to the strong hybridization with the Ni 
atoms, graphene-covered Ni can be an efficient source of 
spin-polarized electrons J^iii In addition, it was also re- 
ported that graphene passivates Ni(lll) against oxygen 
exposure^^ which makes the graphene/Ni(lll) system 
a promising candidate for applications in carbon-based 



magnetic media and spintronic devices. 

In this study, we employ three different magnetic 
imaging techniques to investigate Ni(lll) and graphene 
on Ni(lll), covering length scales from the nanometer 
regime up to several millimeters. With scanning tunnel- 
ing microscopy (STM) the atomic lattice can be resolved 
and bare Ni(lll) is characterized by scanning tunneling 
spectroscopy (STS) where occupied and unoccupied lo- 
cal density of states (LDOS) are probed. Spin-polarized 
STM (SP-STM)i^ can access single domain waUs and 
their response to an applied magnetic field. To investi- 
gate the entire surface domain structure, however, tech- 
niques with a larger field of view are necessary. We used 
scanning electron microscopy with polarization analysis 
(spin-SEM or SEMPA)^^ for vectorial mapping of the 
surface magnetization and Kerr microscopy for field de- 
pendent imaging^ Interestingly, SP-STM and SEMPA 
can access the Ni magnetism through the graphene layer, 
and we show that it is unchanged compared to bare 
Ni(lll). For a deeper understanding of the two sur- 
faces we performed spin-resolved density functional the- 
ory (DFT) calculations. Evaluating the calculated LDOS 
in the vacuum above the surfaces allows the interpre- 
tation of the STM and STS data. As a result of sur- 
face state quenching, we predict an inversion of the spin- 
polarization above the graphene layer with respect to the 
pristine Ni(lll) surface. 



II. EXPERIMENTAL DETAILS 

The same Ni(lll) single crystal (3 mm x 7 mm width, 
1 mm thickness) was used in all experiments. It was 
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cleaned by repeated cycles of 800 eV Ar+ ion etching at 
room temperature (RT) and annealing at T = 1100 K. 
The sample was considered clean when no impurities such 
as carbon or sulfur were detected by Auger electron spec- 
troscopy (AES) and a hexagonal (1 x 1) pattern was ob- 
served by low-energy electron diffraction (LEED). The 
graphene layer was grown on Ni(lll) by chemical vapor 
deposition (CVD)iii Ni(lll) was heated to T = 950 K 
in an ethylene atmosphere (C2H4) of p = 5 x 10^^ mbar 
for 200 s (100 L) and subsequently allowed to cool in 
ultra- high vacuum (UHV). 

The STM experiments were carried out in a multi- 
chamber UHV system with separate chambers for ion 
etching, CVD graphene growth, and STM measurements. 
One microscope operates at RT— while another is ther- 
mally connected to a liquid-He bath cryostat reaching 
T = 8.0 ± 0.5 K and features an out-of-plane magnetic 
field of up to B = ±2.5 TJ^ It is equipped with a 
tip exchange mechanism and for spin-averaged measure- 
ments we used W tips, that were hashed in vacuo to 
approx. T = 2300 K. For our SP-STM studies these tips 
were then coated with about 50 atomic layers (AL) of Cr, 
and annealed at T « 500 K for 5 min. An antiferromag- 
netic tip coating is chosen to minimize the magnetostatic 
interaction between the probe and the magnetic sam- 
ple. Constant-current (/) images (topography) and maps 
of differential conductance (d//dt/) ^^'^^ were measured 
simultaneously with closed feedback loop using lock-in 
technique by adding a modulation voltage JJmod = 25 mV 
to the sample bias U . Single dl /AU spectra were taken 
at specific positions with the tip-sample distance stabi- 
lized at [/stab and /stab before switching off the feedback 
loop. 

SEMPA was used for simultaneous vectorial mapping 
of both orthogonal in-plane magnetization components 
at the surfaceji^ Magnetization images of 20 nm lateral 
resolution were acquired at RT using a primary beam of 
6 nA at 8 keV. The SEMPA instrument is located in a 
separate UHV- facility at 5 x 10~^^ mbar base pressure, 
into which graphene-coated Ni(lll) samples were trans- 
ferred through air. To remove the graphene layer, Ar^ 
ion etching at 600 eV without post annealing was used. 
For contrast enhancement 4 AL of Fe were deposited from 
an e-beam evaporator at 0.2 AL/min. 

We used a full-field Kerr microscope with white-light 
illumination^^ working at ambient conditions to investi- 
gate large areas up to several millimeters on the Ni sur- 
face. Kerr microscopy utilizes the magneto-optic Kerr 
effect (MOKE) to visualize the magnetic surface struc- 
ture of an investigated sample. An arbitrary magneti- 
zation component can be selected for imaging by using 
appropriate apertures in the back focal plane. A lateral 
resolution of 300 nm was achieved and imaging in exter- 
nal fields was possible. 



III. STM RESULTS 
A. Ni(lll) 

Figure[T](a) shows a typical surface area of bare Ni(lll) 
including two monatomic steps. Terrace widths vary with 
lateral position in the range of 20 - 200 nm. As seen 
in Fig. [IKb) the atomic lattice can be resolved on the 
terraces. Fig. [TJc) displays a line profile along a closed 
packed row (white line in (b)) and the interatomic spac- 
ing is in agreement with the nearest-neighbor atomic dis- 
tance of Ni (2.49 A). As expected for closed packed sur- 
faces like fcc(lll), the corrugation is comparatively small 
and lies in the range of 3 - 5 pmj^i 

Despite the nice ordering of surface Ni atoms, residual 
contamination is present in the sample, in particular sub- 
surface defects. These defects have very low corrugation 
in constant current images, but are clearly seen in dl /AU 
maps, e.g. at [/ = +1 V (see the inset in Fig. Hfd)). We 
therefore took care to measure AI/AU spectra on defect 
free areas: Fig. [TJd) shows an average of three spec- 
tra measured at positions as indicated in the inset. We 
observe two broad maxima, 500 meV below and about 
400 meV above the Fermi level, Ey, respectively. The 
spectrum agrees well with the experimental results of 
K.-F. Braun and co-workersi^ Our first-principle calcu- 
lations (Section |V| attribute these features to a minority 
spin surface resonance below and the Shockley state of 
both spin character above E-p. 

To investigate the magnetic properties of the Ni(lll) 
surface we use SP-STM: Fig. [2](a) shows a d//dC/ map 
of a sample area exhibiting nearly horizontal steps. The 
darker and brighter regions indicate magnetic domains 
of Ni(lll) and the dashed lines indicate domain walls. 
To prove the magnetic origin of the observed contrast we 
apply an out-of-plane magnetic field of i? = 50 mT and 
as a result both domain walls have moved to the left as 
seen in Fig. [IJb), the right one by about 150 nm. Line 
sections across the walls, as indicated by a white box in 
Fig. [5Ja) are shown in Fig. [D^c). Fitting a standard wall 
profilo^^ (solid line), 

tanh((x)/(u;/2)), (1) 

where a; is a lateral distance, to the experimental data 
(dots) yields wall widths of w = 88 ± 20 nm and w = 
60 ± 16 nm, respectively. The wall width determined for 
the wall in (b) at 4-50 mT is it; = 62 ± 17 nm. This shows 
that while we can move domain walls in the out-of-plane 
field of -|-50 mT the width of the wall in Fig. [2] is not 
altered within the error of the measured width. The rea- 
son why these walls can be moved by the out-of-plane 
field becomes clear from the volume domain structure as 
deduced from SEMPA measurements fSection lVII[) . The 
magnetic contrast, which can be defined as the asym- 
metry of the d//dC/ signals of bright (b) and dark (d) 
areas 
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FIG. 1. (color online) Ni(lll), spin-averaged data: 
(a) Constant-current image {U = —1 V, 7 = 0.5 nA) taken 
at T = 8 K. (b) Atomically resolved image {U — —4 mV, 
1 = 5 nA) taken at RT. The diamond highlights the unit cell, 
(c) Height profile along the line displayed in (b). (d) The 
dl/dU spectrum (?7atab = +1 V, /stab = 2 nA, C/mod = 
80 mV) is averaged over data taken at three different loca- 
tions marked as red crosses in the d//df/ map (inset). 



A{U) = 



dI/dU{U)b-dI/dU{U)d 
d//d[/(t/)b + d//dC/(C/)d ^ 



(2) 



is typically low on the Ni(lll) surface. Evaluating the 
data in Fig. Hyields a value of only A(-200 mV) = 2 %. 

The spin-resolved dl/dU map in Fig. [2ljd) has a width 
of 1.5 and shows three areas of different intensity, 
i.e. magnetic domains, and an atomic step. The asym- 
metry between the highest and lowest signal amounts to 
2 %. The occurrence of several domains with different 
magnetization directions on this length scale indicates 
already an interesting overall magnetic structure, which 
requires an imaging technique with a larger field of view. 



B. Graphene on Ni(lll) 

The graphene layer is commensurate with Ni(lll) due 
to the small lattice mismatch)^ A typical sample of 
graphene on Ni is shown in Fig. ^a.). The flat terraces in- 
dicate perfect single domain graphene formation without 
any visible defects. At higher magnification in Fig. |3l^b) 
a triangular lattice is seen rather than the honeycomb 
structure of the carbon atoms (the unit cell is highlighted 
by the diamond as shown in the inset). This is not 
surprising since neighboring carbon atoms occupy non- 
equivalent sites on the Ni substrate: in the ball model 
in Fig. EJd) the carbon atoms labeled A reside on top of 
the Ni atoms of the first layer, while the carbon atoms 
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FIG. 2. (color online) Ni(lll), spin-polarized data: Magnetic 
dl/dU maps taken at (a) B = T and (b) B = -1-50 mT. 
Dashed lines show the shift of domain walls in the external 
magnetic field, (c) Line sections across domain walls marked 
by white box in (a). Gray circles and black lines represent 
experimental data and fitted profiles, respectively, (d) Mag- 
netic dl /dU map of 1.5 fim width taken at B = T. All data 
measured at [/ = —200 mV and I = 2.5 nA. 



labeled B are at positions of the Ni atoms of the third 
layer (fee hollow site).— 

The white line in Fig. ^h) indicates the position of the 
line profile in Fig. EIc) and the distance between max- 
ima is the same as for the bare Ni(lll) surface shown in 
Fig. [TJb,c), reflecting the nearest- neighbor distance of Ni 
(2.49 A). However, in contrast to Fig. mb,c), fee and hep 
hollow sites are now distinguishable, due to the B type 
atoms on fee positions. In addition, the atomic corruga- 
tion of 10 - 15 pm for graphene on Ni(lll) observed here 
is roughly a factor of 3 larger than that measured on bare 
Ni(lll). Both a triangular lattice structure and an en- 
hanced corrugation seen in STM images are a purely elec- 
tronic effect originating from graphene pz states around 
E-p and the quenching of Ni surface states as will be dis- 
cussed in Section IVl 

To investigate whether we can still measure a magnetic 
signal on graphene-coated Ni(lll) we perform SP-STM 
measurements. Figures HJa) and (b) show dl/dU maps 
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FIG. 3. (color online) Graphene on Ni(lll), spin-averaged 
data: (a) Constant current overview image (U = —1 Y , I — 
0.5 nA), taken at T = 8 ± 1 K, showing fiat terraces and four 
monatomic steps, (b) Zoom-in on a terrace {U — +2.5 mV, 
7 = 5 nA) showing the atomic lattice, (c) Height profile 
along the line depicted in (b). (d) Top and side views of the 
graphene/Ni top- fee structure: Black color indicates carbon, 
blue color Ni atoms. Red dots indicate the positions of the 
empty spheres (see Section |V|. The unit cell is highlighted 
by the diamond. 



measured in (a) forward and (b) backward scan direc- 
tion with the same bias voltage U = —200 mV as in 
Fig. [2Ia,b). Wc again see areas of different dl/dU con- 
trast as a result of different magnetic domains in the 
image area. This means that the electron density a few 
Angstrom above the surface (at the position of the tip) 
is substantially spin-polarized, despite the fact that the 
carbon atoms are expected to carry a very small mag- 
netic moment.— >2^i^ We conclude that we probe the 
magnetic structure of Ni(lll) under the graphene layer. 
Evaluating the magnetic signal strengths for the dl/dU 
maps in Fig. |4l[a,b), we find a magnetic asymmetry of 
200mV) = 4 %. We observe two types of boundaries 
between homogeneously magnetized areas: while the left 
one can clearly be identified as a domain wall with a 
width of w = 51 ± 6 nm (see line profile and fi t^^'^'' in 
Fig. mjc)), the right one displays a non-continuous tran- 
sition which appears at different lateral positions in for- 
ward and backward scan direction. This sharp transition 
is therefore not a domain wall, but instead it is an ar- 
tifact resulting from magnetostatic interaction with the 
tip (see sketches in Fig. |l{a,b)). Its lateral position de- 
pends on the scan direction (i.e. from left to right or vice 
versa) since the domain wall is pushed along the scan di- 
rection by the tip until it snaps back. This means that 
the tip used in this experiment exhibits a non-negligible 
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FIG. 4. (color online) Graphene on Ni(lll), spin-polarized 
data: Magnetic dl/dU maps measured in (a) forward and 
(b) backward scan directions {U = —200 mV, / = 0.5 nA). 
While the left wall is static, the right one is moved by the 
stray field of the tip (see text), (c) Line section as indicated 
in (a) and a fit of a general wall profile,— yielding a wall 
width of about 50 nm. Data taken at B = T. 



stray field, most likely due to picking up magnetic mate- 
rial from the sample. Such an influence is frequently ob- 
served when using ferromagnetic tips to investigate soft 
magnetic materials both in SP-STM and MFMJ^ Rea- 
sons for the two walls in the image area responding non- 
equivalently might be different magnetization directions 
or different wall types. The fact that the magnetic struc- 
ture of Ni(lll) can easily be changed by small external 
magnetic fields (cf. Fig. [2jb)) suggests that already small 
amounts of ferromagnetic material at the tip apex are suf- 
ficient to observe interactions with the domain structure. 



IV. COMPUTATIONAL METHOD 

To gain a detailed understanding of the observations 
made in our STM experiments we performed DFT calcu- 
lations of the pure Ni(lll) and the graphene/Ni(lll) sys- 
tem using the projector augmented wave^ based Vienna- 
ab-initio-Simulation-Package (VASP)i ^^'^^ We employed 
the generalized gradient approximation (GGA)^ to the 
exchange correlation potential and accounted for van der 
Waals interactions in the calculations involving graphene 
in the framework of the DFT-D2 methodi^ We calcu- 
lated the electronic properties of a clean Ni(lll) slab 
and a Ni(lll) slab coated with graphene on one side both 
containing 15 atomic layers of Ni and ~ 18 A of vacuum 
between periodic images of the slabs. 

The triangular unit cell of the fee Ni(lll) surface and 
the triangular graphene unit cell exhibit a lattice mis- 
match of 1.8 %. In our calculations, we used a unit 
cell with the experimental lattice constant of the Ni(lll) 
surface a ~ 2.49 A^^ and the graphene layer deposited 
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in the so-called top-fcc arrangement (Fig. |3Kd)), which 
has been established as the energetically most favorable 
one in experimental and theoretical invcstigationsj^^^ 
All geometries considered in our calculations were op- 
timized until all forces were below 0.01 eV A-\ For 
obtaining the LDOS the Brillouin zone integrations were 
performed with the tetrahedron method^ using 36 x 36 
k-point meshes. 

We simulated high-resolution STM images by calculat- 
ing the position dependent LDOS in the vacuum and in- 
tegrating the energy window from -100 meV to 100 meV. 
To simulate STM spectra, we have calculated the LDOS 
inside so-called empty spheres placed at 3.6 A above the 
surfaces in the vacuum region of the slab. The LDOS in- 
side the empty spheres was calculated assuming an STM 
tip with an s-wavc symmetric apex state and we care- 
fully checked that the LDOS at lower/higher distances 
(between 2 A and 7 A) from the slab shows a smooth 
trend and yields qualitatively the same values. We con- 
sidered spheres at four different lateral positions above 
the surfaces (see red dots in Fig.[3jd)) to investigate lat- 
eral modulations in the STM spectra. 

In STM experiments, there is an electric field between 
a tip and a sample, which modifies the shape of the tun- 
neling barrier The exact shape of the tunneling barrier 
is unknown but can be approximated by a trapezoidal 
barrier in the most simple model. Besides any density 
of states effects, this leads to an energy and tip-height 
dependence of the d//d[/-signal according to^^ 




where $ is the work function of the tip and the sample 
and s is the distance between the tip and the sample. 
Expanding the exponent to first order in U we arrive at 

-^coexp(^--j, (4) 

where Eq is a constant depending on the materials of the 
tip and the sample as well as on their distance. 

To simulate STM spectra, we therefore use the vacuum 
LDOS from our DFT calculations (which accounts for 
the tunneling barrier due to the sample work function) 
and scale this vacuum LDOS by a factor of exp{—E/Eo) 
to account for electric field induced dependencies of the 
effective tunneling barrier height on the bias voltage. Eq 
is treated fitting parameter. In Section |V] we use 
Ef) = 2 eV, which leads to good agreement of our first- 
principles calculations with the experimental data. 



V. THEORY RESULTS 

We start with the simulated STM images (Fig. [5]) and 
compare them to the experiments (Figs. [IJb) and[3l^b)). 




FIG. 5. (color online) Simulated STM images calcu- 
lated at a height of 3.6 A over (a) clean Ni(lll) and 
(b) graphene/Ni(lll). The orientation of the graphene lat- 
tice is the same as in the unit cell depicted in Fig. [Sj^d): the 
C atoms visible in the image belong to the B sublattice. 



Over the clean Ni(lll) surface, we find a triangular lat- 
tice of protrusions visible as bright spots (Fig. [5ja)), 
which is in good agreement with the experimental STM 
image (Fig. mb)). Our calculations show that the pro- 
trusions are centered above the atoms of the topmost Ni 
layer. 

Our calculations as well as Refs. [13 and [H yield 
graphene adsorbing at a distance of ^ 2.10 A above 
the Ni surface with the graphene layer being almost flat, 
the height difference between the carbon sublattices (a 
structural corrugation) is about 0.5 pm. However, due 
to an electronic effect, the simulated STM images of 
graphene/Ni(lll) also show a triangular lattice struc- 
ture (Fig.[5jb)) again in good agreement with the exper- 
iment (Fig. EJb)). In the convention of Fig. ^d), our 
calculations yield the highest vacuum LDOS above car- 
bon atoms of the B sublattice, i.e. those carbon atoms 
not located above a Ni atom from the topmost Ni(lll) 
layer. This feature is stable with the simulated tip height 
and bias voltage. We thus conclude that the highest pro- 
trusions in the STM images of graphene on Ni(lll) cor- 
respond to carbon atoms in sublattice B. 

In the STM experiments of graphene on Ni(lll), the 
graphene sublattices A and B exhibit an apparent height 
difference on the order of 10 pm. Since the structural 
corrugation is only about 0.5 pm, this enhanced corru- 
gation in the experiments must have an electronic but 
not a structural origin. The electronic states responsible 
for this corrugation will be discussed together with the 
calculated STM spectra below. 

In the experiment we obtained energy resolved 
STM spectra and spin-resolved differential conductance 
(dl/dU) maps, which we now compare to spectra from 
our DFT calculations. The calculated STS spectrum for 
the pristine Ni(lll) surface is shown in Fig. [SJ^a). It ex- 
hibits broad maxima in the energy range between -1 eV 
and -0.5 eV as well as above Ep at energies > 0.5 eV. We 
find that the Ni spectrum does not change with the lat- 
eral position of the tip. A comparison to Fig. [IJd) shows 
that the calculated spectrum is in qualitative agreement 
with the experimental results. 
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FIG. 6. (color online) (a) Simulated STS spectrum at 
3.6 A above the clean Ni(lll) surface, (b) Spin contrast 
Ap (see text) at 3.6 A over the pristine Ni(lll) surface and 
the graphene-coated Ni(lll) surface respectively. 



Our calculations show that the contribution to the 
STS signal arising from majority (spin-up) and minor- 
ity spin states (spin-down) differs clearly, which is well 
in line with spin contrast being achievable in our SP- 
STM experiments. To simulate the magnetic contrast 
provided by the SP-STM we calculated the spin contrast 
Ap = {p^ — p^)/(p-|- -I- Pi) from the spin polarized den- 
sity of states, p-|-|. The spin contrast, averaged over the 
four different spheres in the vacuum (Fig. [3Jd)), for the 
clean Ni(lll) surface and the graphene-coated surface is 
shown in Fig.[6ljb). In agreement with the experiment we 
find a pronounced spin contrast over clean Ni as well as 
over the graphene-coated surface. Thus, the spin polar- 
ization in the vacuum above the surface is not suppressed 
by the graphene coating and the ferromagnetic domain 
structure remains visible in SP-STM. Interestingly, the 
sign of the spin polarization in the vacuum LDOS is re- 
versed for graphene-coated Ni as opposed to pristine Ni 
at energies below -0.23 eV as well as above Ep. This spin 
contrast change may be accessed in future experiments, 
where one would have to perform an SP-STM measure- 
ment on a partially graphene-coated and partially clean 
Ni(lll) sample. 

We now aim at identifying which states of the Ni(lll) 
and the graphcne/Ni(lll) systems give major contribu- 
tions to the tunnel current in the STM experiments and 
thus understand the physics behind the calculated and 
observed STM images, spectra and spin contrasts. To 
this end, we calculated the band structure and analyzed 
the corresponding wave functions of the clean and the 
graphene-coated Ni surface. The band structure of a 
clean Ni slab is shown in Fig. [7] for majority (a) and 
minority (b) spin states. In this figure we use a so-called 
"fat band analysis"— where the displayed thickness of a 
band represents the strength of the property of interest. 
Here the wave function, |^'„,fc), belonging to a given band 
n at a given fc-point is projected onto an s orbital, \L) , 
localized inside an empty sphere at 3.6 A above the re- 
spective surface. The overlap |(^'„^fe|L)p is then depicted 
as the thickness of the corresponding band. 

For clean Ni(lll), the dominant contributions to the 



vacuum LDOS (blue) above E-p arise from upward dis- 
persing bands, which have their energy minima at the 
r point at energies of 10 meV (majority spin electrons) 
and 140 meV (minority spin electrons). These states can 
be characterized as surface states or surface resonances. 
The upward dispersing feature with a minimum at the F 
point for majority spin is the well-known Shocklcy sur- 
face state with mixed Ni and d^z^-r^ character at F 
and dxz.yz admixtures away from FiiSiii A similar fea- 
ture for the minority spin was identified as a resonance 
of mixed Pz and dxz,yz character 

The downward dispersing feature starting at -0.6 cV 
below Ep for the majority states and dispersing along 
the F — >■ K and F — >■ M directions in the Brillouin zone 
was identified as a surface resonance This specific res- 
onance derives mainly from Ni dxz,yz states with a small 
admixture of Ni pz states away from the F point. At 
the F point itself, however, it shows no spectral weight. 
For the minority spin component a similar state begins 
slightly above Ep and disperses down below -1 eV along 
the V ^ K and F ^ M directions. This feature, show- 
ing pronounced weight in the vacuum, was identified as 
a surface resonance with the same orbital characteristics 
as the Shockley state described above4i^ This state is 
responsible for the minority spin polarization dominating 
in the energy region from ^ — 1 eV up to 0.1 eV seen in 
the spin contrast in Fig. [5Ib). The STS spectrum shown 
in Fig. [6ja) also becomes clear now: The features of the 
spin-resolved spectra can be attributed to surface elec- 
tronic features of Ni. The broad peak below Ep stems 
mostly from the downward dispersing feature of the mi- 
nority spin states in Fig. [7l[b), whereas the unoccupied 
spectrum is dominated by the Shocklcy state of both spin 
character. 

Over the graphene-coated surface the situation changes 
qualitatively. Fig. [7l[c,d) show the fat band analysis for 
the Ni slab coated on one side with graphene. Over 
graphene, we find dominant contributions to the spec- 
tral weight in the vacuum (red fat bands) arising from 
free electron like states as can be seen along T K and 
F M directions. The graphene p^-derived bands be- 
tween the K and M points show some smaller weight in 
the vacuum as well. The pronounced surface states and 
resonances seen over the Ni surface cannot be seen over 
graphene, even at very low heights over the surface. The 
band derived from the Shockley state vanishes explicitly 
over the graphene-coated slab surface. Only a band from 
the uncoated surface of the slab remains. Thus, graphene 
quenches the surface resonances and surface states. Since 
the surface states are mainly responsible for the sign of 
the spin contrast in the vacuum over the Ni surface the 
quenching of these states leads to the predicted reversal 
of the sign of the contrast (cf. Fig. [ni[b)). 

We note that this reversal of the spin contrast in the 
vacuum LDOS above graphene does not mean that the 
Ni magnetization is reversed beneath graphene. We find 
in agreement with earlier studies^^ that the magnetic mo- 
ment in the Ni interface layer to graphene is about 20 % 
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FIG. 7. (color online) Band structure of a clean 15-layer slab of Ni(lll) for majority (a) and minority spin (b) components. The 
same slab, coated with graphene on one side, is shown in (c) for the majority and in (d) for the minority case. The projections 
of the bands onto empty spheres in the vacuum at 3.6 A, averaged over their lateral position (red dots in Fig. |3][d)), are shown 
as fat bands. Here, the vacuum projections above the Ni are indicated by blue color, over graphene by red color. In (e) and 
(f) the contributions of each band to the corrugations measured in STM are visualized as fat bands (see Eq. ((5}) in the same 
color code. For visualization purpose, the corrugation above Ni has been enhanced by a factor of 4. 



smaller (0.51 /ie) compared to the bulk value of 0.65 /iB. 
Additionally, a small spin magnetic moment of -0.02 /ie 
for Ca and 0.03 /xb for Cb atoms respectively is induced 
in the graphene layer 

We finally address the question why the apparent 
height corrugations measured in STM (see Section 
are larger in the graphene/Ni system than on the bare 
Ni(lll) surface. To this end, we compared the vacuum 
amplitudes of the states of clean Ni(lll) and graphene- 
coated Ni systems at different lateral positions within 
the unit cell (red dots in Fig. |31Id)). The variation of 
the state amplitude with the lateral position is visual- 
ized in Fig. [He,f): for each band |^'Ti,fc}, the thickness 
d of the colored curves (blue: Ni, red: graphene) corre- 
sponds to the standard deviation of the vacuum projec- 
tions |(^'r,.fe|ir)P with lateral sphere position r: 



The line thicknesses are normalized to the (Brillouin zone 
integrated) local density of states p^i{E) at the respective 
energies. 

The variation of the wave function amplitudes with 
the lateral sphere position is clearly higher over graphene 
than above clean Ni, which can be seen in Fig.[7|e,f). It is 
visible that the corrugations measured above graphene on 
Ni mainly arise from graphene derived states. Similar 



to the case of graphite)^ the peculiar symmetry of these 
states induces an asymmetry in STM images. For small 
bias voltages, hence close to Ep, the upper graphene Pz- 
derived band crossing Ep near the M point contributes 
strongly to the corrugation in the majority spin chan- 
nel. The lower pz band and small contributions from 
Ni d bands induce the corrugation in the minority spin 
case. Above clean Ni(lll) our STM experiments measure 
smaller corrugations, which mainly originate from slight 
lateral variations in the surface states near the center of 
the Brillouin zone. It becomes clear that the corrugations 
above the graphene sheet are of electronic origin and are 
also brought about by the quenching of the Ni surface 
states as can be seen in our data. 



VI. KERR-MICROSCOPY RESULTS 

To identify the large-scale magnetic domain structure 
we studied the domain pattern of Ni(lll) by means of 
Kerr microscopy. An in-plane measurement shows a re- 
curring magnetic pattern that varies slightly for different 
demagnetization cycles while we found no indication for 
an out-of-plane component. In Fig. |5] the horizontal (a) 
as well as the perpendicular component (b) of the in- 
plane magnetization is shown. The in-plane magnetic 
pattern is characterized by magnetic structures of two 
length scales: we observe a stripe domain pattern with a 
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FIG. 8. Kerr microscopy images of the magnetic structure of 
the Ni(lll) single crystaL The gray-scale amplitude of the 
images is proportional to the in-plane magnetization compo- 
nent which is indicated by a double arrow in the lower right 
corner. 

stripe width in the range of 3 /im to 6 fxin which varies 
the orientation in different sample areas. Inside each of 
these stripes one can see a wavy pattern indicative of a 
magnetic fine structure on a smaller length scale. A real- 
time observation of the Ni(lll) single crystal during the 
application of an external magnetic in-planc field shows 
that at about 13 mT the domain walls begin to move. 
A field of 140 mT is sufficient to create a single domain 
state. 



VII. SEMPA RESULTS 

To get a more detailed picture of the surface magnetic 
domain structure we performed high-resolution SEMPA 
measurements. The SEMPA image in Fig. IH^a) shows 
a characteristic section of the magnetization pattern of 
the graphenc-covcred Ni(lll) single-crystal surface, af- 
ter transfer under ambient conditions. Although no fur- 
ther treatment of the surface has been performed, we ob- 
serve a meaningful magnetic contrast, i.e. a polarization 
asymmetry of 1.2 %. The observation of magnetic con- 
trast in SEMPA without in situ cleaning of the sample is 
most unusual, because of the sub-nm surface sensitivity 
of secondary-electron spin polarization. 

Therefore we must deduce that the Ni(lll) surface is 
effectively passivated by the graphene layer, which is in 
accordance with recent spectroscopy dataJ^ To check for 
any graphene-induced change of signal and/or domain 
pattern Figure IHJb) shows the same area of the sample 
as in (a) after argon ion sputtering. While the domain 
structure is unchanged, the clean Ni surface exhibits a 
much stronger magnetic contrast corresponding to 2.4 % 
asymmetry. The sign of the secondary-electron spin po- 
larization is preserved, which indicates, that in the SE 
cascade process the graphene does not cause a polariza- 
tion inversion. This should not be confused with the cal- 
culation given in Fig. [BJ as SEMPA detects free electrons 
in vacuum, which have an energy > 5 cV with respect 
to the Fermi level. The image quality can be improved 



by depositing a small amount of iron 4 AL) onto the 
surface. Due to its higher saturation magnetization and 
the sub-nm surface sensitivity of SEMPA the iron acts as 
polarizer for the secondary electrons. The result can be 
seen in Fig. |9jc) and again the domain pattern is identi- 
cal to (a) and (b), while the contrast is enhanced to give 
a polarization asymmetry of 6.7 %. These differences 
of the magnetic contrast and thus signal-to-noise ratio 
arc highlighted by the gray-scale images in the bottom 
right corner, which display the experimental results on 
the same scale. 

The SEMPA measurements confirm the results from 
Kerr microscopy and provide images of the surface mag- 
netic domain pattern at higher resolution: it consists of a 
larger length-scale stripe pattern with a width from 3 /im 
to 6 /im. However, these larger stripes are not single do- 
mains but instead the magnetization within the stripes 
is more or less regularly modulated by a second type of 
stripes, on a smaller length-scale from 1 /xm to 3 /im. 
From one of these small stripes to the next, the mag- 
netization changes by approximately 60 degrees, which 
results in a net magnetization of the large stripes. The 
transitions between the smaller stripes are very broad 
yielding a rather wavy pattern without sharp domain 
walls. In contrast, the large stripes are separated by 
sharp Neel-like walls, which are of the 180° type. Fig- 
ure llOf a) shows a higher magnification of the domain 
structure around such a head-to-hcad domain wall. The 
y-component of the magnetization is shown in (b) . The 
striped line indicates the position of the wall profile plot- 
ted in Fig.lTUfd). It can be described by Eq. (11} , and a fit 
yields a domain wall width w = 330 ± 6 nm. The obser- 
vation of much narrower domain walls by SP-STM does 
not contradict the SEMPA results: in first order approxi- 

— 1/2 

mation the width of a Neel wall is proportional to ' , 
where Ki is the first order ma gnet o crystalline anisotropy 
constant. From Refs. and 113 it is known that Ki is 
strongly temperature dependent between 300 K and 4 K: 
-0.0045 MJ/m^ at RT compared to -0.12 MJ/m^ at 4 K. 
Therefore, a reduction of wall width with temperature is 
expected. Using the width from the SEMPA measure- 
ments as starting point, we estimate a wall width at 4 K 
of approximately 70 nm, in reasonable agreement with 
the waU width found using SP-STM. 

To reveal the origin of the observed surface magnetic 
pattern it is crucial to understand why the Neel walls 
have a head-to-head or tail-to-tail configuration as seen 
in the SEMPA images. If this pattern persisted into the 
volume, it would imply a huge amount of dipolar energy. 
Instead, this pattern is indicative of a so-called V-line 
structure^ where two volume domain walls with differ- 
ent orientations merge at the surface into a single line. 
In a cross-section perpendicular to the line, this struc- 
ture appears like a V where the magnetization of the 
center domain collects all the flux that originates from 
the oppositely magnetized side domains (sec Fig. [TUT c)). 
So the surface pattern is flux compensated in the vol- 
ume to reduce the dipolar energy. As the magnetic sur- 
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FIG. 9. (color online) SEMPA images of the magnetic structure of the Ni(lll) surface at the same position, but with different 
surface preparations. Encoded in color is the direction of the in-plane magnetization, as defined by the 360° color wheel. In 
(a) the surface is covered with graphene. (b) gives the signal from the clean Ni surface following sputter-cleaning. For contrast 
enhancement, in (c) a thin iron layer has been deposited. To emphasize the changes in contrast, the lower right part of each 
image gives the y-component of the magnetization in a gray-scale representation. 



face structure of the V-lines imaged in Fig. [10] shows 
only an in-plane magnetization we can interpret it as the 
Neel cap of a V-line, in analogy to the well-understood 
Neel cap of a Bloch wall e.g. in Fe(001)4^ Indeed, our 
measurements give no indication of an out-of-plane com- 
ponent of the magnetization in the domains or the do- 
main walls, in agreement with the Kerr-microscopy mea- 
surement. Thus we conclude that the magnetization at 
the surface is entirely in-plane. This may be surprising, 
as none of the magnetically easy (lll)-directions of the 
Ni crystal lies within the (111) surface of the sample. 
This finding can only be explained as a consequence of 
a reduction of stray field energy at the expense of local 
anisotropy within a certain depth of subsurface volume. 
Estimating the upper limit of the shape anisotropy of the 
crystal with and comparing this to the magneto- 

crystalline anisotropy constant in first order Ki at RT we 
obtain a ratio of 34:1. This might explain the in-plane 
orientation of the observed magnetic pattern. 

The above-described features of the magnetic pattern 
of the surface of the Ni(lll) crystal can be understood 
qualitatively using the quasi-domain branching approach 
for large crystals with strongly misoriented surfaces given 
by Hubert and Schaferi^ For a Ni platelet with (111) sur- 
faces the model assumes 180°-oriented base domains in 
the volume, which reduce the magnetostrictive energy of 
the crystal. The lateral extent of these domains would 
be responsible for the length scale of the larger stripes 
we find on the surface. In the mentioned quasi-domain 
branching concept, quasi-domains with a net magnetiza- 
tion parallel to the surface are introduced, which close the 
flux of the basis volume domains. They are composed of 
alternating domains oriented along the easy directions. 
Each of these first-level branching domains acts as ba- 
sis domain for second-level branching and thus forms its 
own closure quasi-domain at the surface to further re- 



duce the stray-field energy. In the model, the energy 
gain by branching of the closure domains in comparison 
to the amount of domain wall energy needed for further 
branching determines the branching depth that is finally 
observed. In our measurements, we observe a two-level 
branching only, where at least the surface of the second 
level domains is already fully in-plane oriented. This is 
in contrast to the model of Hubert, where quasi domains 
along the out-of-plane canted easy axes are expected even 
on the final level of branching. 



VIII. SUMMARY 

Finally, we would like to summarize the main as- 
pects of this study, (i) We probed the differential con- 
ductance of bare Ni(lll) by STM and in comparison 
to our DFT calculations attributed the observed fea- 
tures to a minority spin surface resonance below and the 
Shockley state of both spin character above Ep- (n) For 
graphene/Ni(lll), STM images showed a triangular lat- 
tice of the atomic structure and an enhanced corruga- 
tion compared to pure Ni(lll). DFT showed that both 
properties are a purely electronic effect originating from 
graphene Pz states around Ep and the quenching of Ni 
surface states, (iii) Kerr microscopy and SEMPA mea- 
surements revealed an entirely in-plane magnetic pat- 
tern of Ni(lll) which stems from a two-level branch- 
ing of the domain structure, (iv) We could easily move 
domain walls in low external magnetic fields or by the 
stray field of a magnetic tip. A single in-plane domain 
state was observed at RT in an in-plane magnetic field 
of 140 niT. (v) The reactive properties of Ni(lll) are 
passivated by the graphene layer well enough for mag- 
netic imaging after transfer through air. The magnetic 
structure of graphenc-coatcd Ni is unchanged compared 
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to that of bare Ni(lll). (vi) As a result of surface state 
quenching, our DFT calculations predict an inversion of 
spin-polarization above the graphene layer with respect 
to the pristine Ni(lll) surface. 
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FIG. 10. (color online) High resolution SEMPA images of 
the Neel-like walls at the Ni(lll) surface. In (a) the in-plane 
distribution of the magnetization of the Ni(lll) surface is 
shown. Encoded in color is the direction of the in-plane mag- 
netization, as defined by the 360° color wheel. In addition the 
direction is illustrated by the arrows. In (b) the y-component 
of the magnetization is shown in a gray-scale representation. 
The striped line is indicating the path of the wall profile shown 
in (d). (c) Sketch of a section across a V-line on the surface, 
illustrating the corresponding volume domain structure. The 
striped arrows indicate the branched structure of the corre- 
sponding domains, (d) The wall-profile is fitted by Eq. U]) 
This fit yields a wall width w = 330 ± 6 nm. 
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